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We analyze the atom field-effect transistor scheme [J. A. Stickney, D. Z. Anderson and A. A. 

Zozulya, Phys. Rev. A 75, 013608 (2007)] using the standard tools of nonequlilibrium dynamics. 

In particular, we study the deviations from the Eigenstate Thermalization Hypothesis, quantum 
fluctuations, and the density of states, both ab initio and using their mean-field analogues. Having 
fully established the quantum vs. mean-field correspondence for this system, we attempt, using a 
mean-field model, to interpret the off-on threshold in our transistor as the onset of ergodicity—a 
point where the system becomes able to visit the thermal values of the former integrals of motion 
in principle, albeit not being fully thermalized yet. 


I. INTRODUCTION 

Similarly to the conventional electronic transistors, in 
atomtronics, a transistor is a nonlinear device where a 
small atom current or atom number controls a large cur¬ 
rent. The existing proposals for an atom transistor can be 
divided into open and closed architectures. The former 
m is conceptually close to the conventional electronic 
devices (the work [ 2 ] also contains a design for a diode). 
The latter, in turn, can be separated into the schemes 
based on an adiabatic population transfer protocol ms] 
(along with the similar atom diode proposals [3 [5]), and 
the schemes where the (large) base current is induced by 
the difference in chemical potentials between the “elec¬ 
trodes”. The base current can be controlled by either 
the internal state of a localized impurity [SHII] or by 
the number of (strongly interacting) atoms in the middle 
site [nna. This latter scheme is nothing else but an 
atomic version of the field-effect transistor (FET). Here, 
the large atomic current from the “source” to the “drain” 
is controlled by a a small number of atoms in the “gate.” 
The interatomic repulsion in the gate is so strong that 
the chemical potential is large even when the number of 
atoms in the gate is small. 

Here we show that the off-on threshold in our transistor 
corresponds to the onset of ergodicity—the point where 
the system becomes able to visit the thermal values of the 
former integrals of motion in principle, albeit not being 
fully thermalized yet. 

II. SYSTEM OF INTEREST 

The paradigmatic atom transistor is realized in a tight 
binding model for three coupled bosonic wells on a line, 
the “source”, the “gate”, and the “drain” [T^- The 
Hamiltonian of the system reads 



where the indices “s”, “g”, and “d” stand for source, gate 
and drain, respectively, e,, and are the site-dependent 
one-body energy and the strength of the on-site two-body 
interactions for a site ry = s, g, d, and Jrjri' is the hopping 
constant between sites 77 and ry' . 

In general, the purpose of the atom FET transistor 
is to control the current of atoms from the source to 
the drain using small variations in the chemical poten¬ 
tial of the gate, pig = Sg — UgUg, controlled in turn by 
its population, Ug. According to the original proposal 
the transition of atoms between the source and the 
gate plays the role of a “bottleneck,” activated only when 
the source and gate chemical potentials become close; in 
contrast, the gate-drain link is made insensitive to the 
gate-drain chemical potential difference, transmitting ev¬ 
ery atom that happens to appear at the gate. This can 
be achieved by providing a comparatively large on-site 
interaction strength for the gate, subsequently detuned 
from the energy of the source. One then chooses the 
hopping constants and the on-site energies in such a way 
that the source-gate transition exhibits a narrow reso¬ 
nance, and the gate-drain transition a broad one. We 
will be using the following representative set of param¬ 
eters: Jsg = —0.1, UgN = 100, Eg = —1.3, Jgd = —1.0, 
Ed = 0.5 and N = 30, with the rest of the parameters set 
to zero: UgN = Eg = U^N = 0. 

Figure [TJi shows the energy spectrum Ea of the tran¬ 
sistor for the chosen set of parameters. Here and below, 
a is the index of the eigenstate |a) of the Hamiltonian 
in Eq. Q: H\a) = Ea\a). The spectrum, bounded from 
both below and above, contains 496 eigenstates. Note 
that the density of states decreases with energy. On the 
upper end of the spectrum, the dominant contribution 
to the energy is provided by the interactions between the 
atoms in the gate, ng; this energy increases quadratically 
with rig, leading to the decreasing (with energy) energy 
spacings. 

Most of the spectrum turns out to be within the limits 
of applicability of the semiclassical approximation. To 
verify this assertion, we check that Weyl’s law holds in 
our system. According to Weyl’s law, if one computes 
how the number of quantum states below an energy E 
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FIG. 1: (color online) Energy spectrum of the atom tran¬ 
sistor. Parameters are Jgg = —0.1, U^N = 100., £g = 
— 1.3, Jgd = —1., £d = 0.5, while U^N = Ss = U^N = 0. 
(a) The ab initio energy spectrum (496 eigenstates) for 
an atom transistor with = 30 atoms. (b) The mean- 
field emulation of the quantum energy spectrum. 496 
initial conditions were uniformly distributed through the 
whole available phase space. The vertical axis gives the 
temporal average of total energy, for each realization. 
The realizations were permuted in such a way that the 
energies increase along the list. The vertical axis reflects 
the position of the realizations in this list. 

depends on the energy E, and then takes the smooth 
envelope of this dependence, the result is proportional to 
the classical phase space volume occupied by these states: 

(# of states with Ea < E) 
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(f’v ~ 0(a:) is the Heaviside step- 

function. Figure shows the inverse of the depen¬ 
dence computed using a Monte-Carlo method: 496 
phase-space points were distributed uniformly within 
an intersection of the whole available volume of the 
phase space and a narrow shell corresponding a con¬ 
stant window of the values of the norm (where norm = 

+ 4 + 

4)), around norm = N. Notice that for this calculation, 
we have chosen as many Monte-Carlo realization as there 
are quantum eigenstates of the system. A priori, one 
number is not related to the other, and more realizations 
would produce a smoother semi-clalssical curve. How¬ 
ever, it can be argued that Fig. [^3 can serve as a semi- 
calssical emulation of the quantum spectrum Fig. 
where the fluctuations of the eigenenergies around a 
smooth envelope have at least the same order of mag¬ 
nitude as their semiclassical counterparts. Indeed, on 
the one hand, the Monte-Carlo method produces energy 
values through a Poisson process; on the other hand, it 
is known that in quantum integrable systems, eigenen¬ 
ergies are distributed as if they are the outcome of a 
Poisson process [Til US]. In contrast, quantum-ergodic 
systems with time-reversal invariance exhibit the type-1 
Wigner-Dyson statistics [IlIISl- In that case, the stan¬ 
dard deviation for the spacing between two neighboring 
levels is l/\/(4/7r) — 1 « 2 times lower than in the case 
of Poisson statistics. 


III. ATOM TRANSISTOR FROM THE POINT 
OF VIEW OF QUANTUM NONEQUILIBRIUM 

DYNAMICS: THE OFF-ON CROSSOVER AS 
THE POINT AT WHICH THE SYSTEM 
BECOMES ABLE TO REACH THERMAL 
EQUILIBRIUM VALUES 

Figure]^ explores the central object in the eigenstate 
thermalization theory [TMT^ —the expectation values of 
the relevant observables; in our case, this is the relative 
occupation of the drain: (a|hd|Q:)- Notice that this ex¬ 
pectation value, when plotted as a function of a, does 
not collapse to a single line; this means our system has 
not reached eigenstate thermalization. Let = Njj/N, 
where = {b\^b^). Figure]^ shows the classical 

infinite time averages of the drain occupation, 

{nd)t= lim -- / dtudit) , 

imax—Cniax Jo 

where the pseudo-eigenstates of Fig. Et are used as the 
initial conditions. The overall behavior of the quantum 
and classical expectation values stand in good correspon¬ 
dence (see comment 24 in Ref. |2()jf. 

Quantum fluctuations of Ud in the eigenstates ja). 


is the corresponding mean-field Hamiltonian, expressed a[hd]a = \/{a\{hd)'^\a) - {{a\hd\a))‘^ , 

through the action-angle variables of the zero-hopping 

analogue of the Hamiltonian in Eq. Q, 4 = shown in Fig. In accordance with the conjecture 
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FIG. 2: (color online) Expectation values of the drain oc¬ 
cupation. (a) Quantum expectation values in the eigen¬ 
states. The solid horizontal line reflects the average over 
all the eigenstates. The system parameters are the same 
as in Fig.j^a). (b) Time average of the drain occupation 
of the atom transistor, in the mean-field approximation, 
for all the realization used. The long solid horizontal 
line reflects the average over all the realizations. The 
short solid horizontal line is the average over 25 real¬ 
izations with the lowest energies; these occupy approx¬ 
imately 1/20 of the total available phase space volume. 
This portion of the phase space constitutes the domain 
of operation of the atom transistor (see Figs. |3and|^. 
The system parameters are the same as Fig.fiTb) 


FIG. 3: (color online) Fluctuations of the drain occupa¬ 
tion. (a) Quantum fluctuation in the eigenstates. The 
solid horizontal line reflects the average over all the eigen¬ 
states. The parameters are the same as in Fig. [^a). (b) 
Temporal variance for the initial states generated, in the 
mean-field approximation. The solid horizontal line re¬ 
flects the average over all realizations. The parameters 
are the same as in Fig. j^b) 

The conventional figure of merit is the current to the 
drain, characterized by the initial slope in the depen¬ 
dence of rig vs. time (in our case, the slope of the red line 
between t — 0 and t « 20, in Fig.[^). Below, we will use 
an integral figure of merit—the infinite time average of 


expressed in Ref. EH, the classical temporal variance, 

O'* [rid] = \l ( lim - ((nd)t)^ , 
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shown in Fig. constitutes a faithful classical analogue 
of the quantum fluctuations. 

Figure]^ shows the detuning of the source-gate transi¬ 
tion as a function of the relative gate occupation. The 
resonance occurs at rig = 0.013. Recall that in our case, 
/is = Es and /ig = Eg — UgUg. 

In all numerical experiments below, the initial value of 
the drain occupation rid is zero or very close to zero. 

At zero rig initially, only a minimal conductance is ex¬ 
pected, and this is what we observe (Fig.[^). To the con¬ 
trary, at the resonant point, ng(t = 0) = 0.013, the tran¬ 
sistor is supposed to provide the maximal source-drain 
conductance, and this is indeed what we find (Fig. 



FIG. 4: (color online) Detuning between the gate and the 
source chemical potentials, as a function of the relative 
occupation of the gate, Ng/N. The vertical half-width of 
the line equals the coupling constant t4g between source 
and the gate. The rest of the parameters is the same as 
St Fig.[|;b). 
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FIG. 5: (color online) Dynamics of the source (solid, 
blue), gate (short-dashed, green) and drain (dashed, red) 
occupations in the “on” and “off” modes, (a) “Off” 
mode (low values of Ug = Ng/N in Figs. and |^. (b) 
“On” mode (value Ug = Ng/N close to a resonant value 
of 0.013, according to Figs. [^and[^. 

the drain occupation, (nd)t—to facilitate the analysis of 
the system from the point of view of ergodicity or its 
absence. 

In Fig.l^ we trace the (nd)t vs. ng(t = 0) dependence, 
along with the temporal variation of n^. The slope of 
this dependence, /3 « 16, we identify as gain. Unexpect¬ 
edly, we find that even at the resonant point, the aver¬ 
age drain occupation does not reach its ensemble average 
(also represented by the short horizontal line in Fig. [^. 
Instead, as the system approaches the resonant point 
Ugit = 0) = 0.013, the ensemble average starts falling 
within the range of the temporal fluctuations; not enough 
to be predominantly within the range, but enough for this 
effect to be detectable. 


IV. CONCLUSION 

In this work, we studied the atom FET transistor 
scheme suggested in Ref. [nma using the tools of quan¬ 
tum non-equilibrium dynamics |20j . We first justified 
the applicability of the semiclassical approximation ap¬ 
plied to the standard measures of the non-equilibrium. 
We then focused, using a semiclalssical model (equiva¬ 
lent to a mean-field one), on the initial conditions with 
zero drain occupation, and used the gate occupation as 
the knob that controls further dynamics. Instead of the 
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FIG. 6: (color online) Time average of the drain occupa¬ 
tion of the atom transistor (dashed line), along with its 
temporal standard deviation (error bars), in the mean- 
field approximation, as a function of the initial occupa¬ 
tion of the gate. 1015 initial conditions were uniformly 
distributed inside a phase space volume corresponding to 
the lowest (in energy) 1/20 of the total available phase 
space volume. The corresponding window of energies is 
between E^i^/N « 0.08 and E^^y^/N « 0.13, out of the 
full available energy range of Uniin,fuii/./V = Ey^iyi/N « 
0.08 and Emax.fuii = 47. Subsequently, only the ini¬ 
tial conditions corresponding to low initial conditions 
{N^{t = 0) < 0.05A^) of drain occupation were selected. 
A linear fit (dashed line) gives an estimate of the gain, 
/? = {Nd)t/Ng{t = 0), as /3 « 16. Solid horizontal line re¬ 
flects the average over all the realizations in the window, 
with no selection of the initial values of the drain occupa¬ 
tion Ad. Note that only when the initial gate occupation 
approaches the resonant value of Ug = .013 (see Fig. 13, 
the mobility in the phase space becomes sufficient for the 
drain occupation to reach its thermal value (solid line), at 
least at some instances of time. The system parameters 
are the same as in Figj^b) 

traditional source-to-drain current, as the transistor out¬ 
put we choose the infinite time average of the drain occu¬ 
pation: this change brings us closer to the conventional 
measures of non-equilibrium. We then studied how the 
drain occupation depends on the initial gate occupation. 
At zero gate occupation (transistor “off”), no atoms are 
transmitted to the source. To the contrary, at a point 
where the gate’s chemical potential levels with that of 
the source (transistor “on”), the drain becomes popu¬ 
lated. One may expect that at this point, the transistor 
is fully thermalized, i.e. that the infinite time averages of 
the observables approach their (microcanonical) ensem¬ 
ble averages. However, it turned out that even at the 
resonant point, the drain occupation is still far below its 
thermal value. Instead, we found that the “on” regime is 
characterized by the accessibly of the thermal values of 
the drain in “principle”. Namely, around the “on” value 
of the initial gate occupation, the temporal fluctuations 
of the drain occupation become capable of “touching” the 
thermal values for some periods of time, without spend¬ 
ing much time there. The overall conclusion is that the 
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ergodicity is not a necessary condition for the operation 
of atom transistor (while the conventional semiconductor 
devices do operate close to a thermal equilibrium): this 
conclusion may allow one to broaden the search for an 
optimal configuration of an atom transistor. 
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